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SUMMARY 


The  objective  of  this  program  on  coherent  optical  adaptive  techniques 
(COAT)  is  to  investigate  algorithms  and  techniques  that  can  reduce  the  beam 
distortions  caused  by  thermal  blooming.  This  report  covers  the  contract 
period  from  1  November  1975  through  31  Mary  1976.  The  work  during  this 
period  has  centered  on  the  construction  and  preliminary  testing  of  a  deform¬ 
able  mirror  and  associated  electronics. 

The  construction  of  a  37-element,  all-beryllium  deformable  mirror 
has  been  completed.  Final  polishing  of  the  mirror  is  still  not  complete 
because  of  the  inability  of  a  supplier  to  deliver  the  bulk  of  the  37  piezoelectric- 
stack  driver  elements.  Only  19  of  the  37  stacks  have  been  installed  at  this 
writing.  Some  preliminary  tests  on  the  mirror  have  determined  its  excur¬ 
sion  sensitivity,  influence  function,  and  frequency  response.  At  present, 
the  mirror  delivers  only  0.28  jxm  of  surface  motion  for  150  V  of  drive, 
slightly  less  than  one  half  the  design  goal.  The  frequency  response  shows 
some  unusual  resonance  behavior,  but  is  usable  to  at  least  30  kHz.  The 

1  -  5 

influence  function  has  a  shape  given  by  IN(r)  =  exp(-Br  )  and  has  a  value 
of  IN  =  0«23  at  the  center  of  the  adjacent  actuator.  The  faceplate  will  be 
cut  thinner  to  reduce  this  value  to  about  0.  15. 

The  computer  simulation  studies  of  irradiance  tailoring  reported  in 
the  previous  report  were  continued  to  correct  an  error  in  the  results  pre¬ 
sented  there.  The  basic  conclusion  is  unchanged,  however:  the  best  way  to 
reduce  thermal  blooming  is  to  increase  the  transmit  aperture  size  and 
make  the  beam  as  uniform  as  possible. 

We  have  completed  construction  of  high-voltage  mirror  drivers  and 
a  dither  injector  for  use  with  the  deformable  mirror.  Construction  has 
begun  on  a  14-function  Zernike  polynomial  generator.  Computer  simulation 
studies  have  identified  which  Zernike  polynomials  make  the  greatest  con¬ 
tribution  to  blooming  compensation.  A  computer  code  has  been  developed 
that  provides  data  on  generating  Zernike-polynomial  mode  shapes  with  a 
finite  number  of  deformable  mirror  elements  with  minimum  mean  square 
error  between  the  actual  Zernike  shape  and  the  mirror-produced  surface. 
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During  the  next  period  the  deformable  mirror  construction  and 
characterization  will  be  completed.  The  Zernike-polynomial  generator 
will  be  completed,  and  performance  tests  will  begin  with  the  deformable 
mirror  in  the  DARPA/RADC  COAT  system* 


Publication  Review 


This  technical  report  has  been  reviewed  and  is  approved. 


Robert  F.  Ogrodnik 
RADC  Project  Engineer 
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I. 


INTRODUCTION 


A.  Program  Objectives 

The  primary  objective  of  this  program  on  coherent  optical  adaptive 
techniques  (COAT)  is  to  analyze  and  experimentally  demonstrate  adaptive 
multidither  correction  algorithms  that  can  reduce  beam  distortions  caused 
by  thermal  blooming.  The  use  of  fixed  transmitter  intensity  profiles  for 
reducing  thermal  blooming  will  also  be  investigated. 

B.  Research  Program  Plan 

The  research  program  utilizes  the  21 -channel  DARPA/RADC  experi¬ 
mental  COAT  system  built  and  tested  on  contracts  F30602-73-C-0248  and 
F30602-75-C-0001.  Computer  simulation  codes  developed  on  these  con¬ 
tracts  and  on  other  programs  (e.g. ,  NSWC  contract  N60921 -74-C-0249) 
have  been  used  for  the  analytical  portions  of  this  contract.  These  codes 
model  the  operation  of  several  types  of  COAT  servomechanisms  as  well  as 
the  time-dependent  propagation  of  optical  beams  in  an  absorbing  and  turbu¬ 
lent  medium.  The  experimental  investigations  require  construction  of  a 
deformable  mirror  for  the  COAT  system  as  part  of  the  program. 

To  accomplish  the  objectives  of  this  contract,  a  12-month  research 
program  consisting  of  three  major  tasks  was  developed.  Task  1  is  an 
analytical  task  and  is  largely  completed.  Task  2  provides  for  the  construc¬ 
tion  of  a  deformable  mirror  and  is  also  nearly  complete.  In  Task  3  we  will 
use  this  mirror  to  study  thermal  blooming  compensation  with  zonal  multi¬ 
dither  control  and  with  modal  control  using  Zernike-polynomial  modes. 

Amendment  No.  1  to  the  contract  has  provided  funds  to  characterize 
the  beryllium  deformable  mirror  being  built  and  to  assess  the  applicability 
of  this  design  to  high-power  laser  COAT  systems.  The  characterization  of 
the  mirror  will  occur  prior  to  the  start  of  Task  3. 
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II. 


TECHNICAL  ACCOMPLISHMENTS 


Ac  Zernike  Polynomial  COAT  (ZEP-COAT) 

At  the  beginning  of  this  program,  we  had  some  hope  that  multidither 
control  of  Zernike  polynomial  functions  (modal  control)  might  produce  a 
greater  degree  of  compensation  for  thermal  blooming  than  localized  (zonal) 
control  of  the  transmitted  phasefront.  As  a  result  of  our  computer  simula¬ 
tion  work  and  experimental  results  at  Hughes  and  at  Lincoln  Laboratory, 
and  through  a  better  understanding  of  the  limitations  of  phase  compensation 
of  thermal  blooming,  we  no  longer  believe  that  modal  control  has  any 
inherent  advantage  over  zonal  control  in  the  degree  of  compensation  that  it 
can  produce.  This  conclusion  is  certainly  true  for  turbulence  compensa¬ 
tion  using  Zernike  polynomial  modal  control,  as  shown  by  the  results  in 
Figure  1.  For  turbulence  compensation,  it  is  the  number  of  control  chan¬ 
nels  (number  of  degrees  of  freedom)  that  counts,  not  what  functions  are 
controlled. 

For  thermal  blooming,  however,  it  may  be  possible  to  obtain  sig¬ 
nificant  blooming  correction  with  far  fewer  control  channels  than  v/ith  a 
zonal  control  system.  The  results  of  Figure  1  also  indicate  that  signifi¬ 
cant  turbulence  correction  is  possible  with  a  limited  number  of  Zernike 
control  channels.  With  fewer  channels  the  complexity  of  the  control  sys¬ 
tem  is  reduced,  and  faster  servo  response  is  possible  when  the  dither 
mirror  response  limits  the  servo  bandwidth.  This  last  conclusion  results 
from  the  computer-determined  result  that  the  maximum  servo  bandwidth, 

(fs)  ,  for  a  given  number  of  sine/cosine  dither  channels,  N  ,  is  given 
max  c 

by 


(fs) 


max 


f 


max 

10  +  1.6  (N 

c 


2) 


(1) 


where  f  is  the  maximum  dither  frequency  allowed  by  the  deformable 
inax 

dither  mirror.  Note  that  N  in  Eq.  (1)  is  the  number  of  control  channels, 

c  - 
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Figure  1.  Analytical  comparison  of  turbulence  compensation  with 
Zernike  polynomials  and  zonal  COAT  control.  The 
transmitter  diameter  is  D^. 
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e 


not  the  number  of  dither  frequencies  (one  half  of  N  in  a  sine/cosine 

c 

system),  and  is  assumed  to  be  a  multiple  of  2. 

In  principle,  a  multidither  COAT  system  can  operate  with  any 
function  set  that  is  orthogonal  (or  nearly  so)  over  the  transmitting  aper¬ 
ture.  That  is,  the  correction  phase,  (j)^,  can  be  generated  according  to 

N 

‘*’c  =  X^(^n  +  ^n  V) 

n=l 

where  the  two-dimensional  functions  f^  are  orthogonal  over  the  correction 

aperture,  a  are  the  correction  amplitudes,  and  6  and  to  are  the  dither 

^  th  ^  ^ 

amplitude  and  frequency  for  the  n  servo  channel.  For  a  zonal  system, 

the  functions  f  are  the  influence  functions  for  the  N  deformable  mirror 
n 

actuators. 

For  optical  systems  with  circular  apertures,  a  logical  choice  of  an 
alternate  orthogonal  correction  set  is  the  Zernike  polynomials.  These  func¬ 
tions  are  orthogonal  over  a  circular  aperture,  and  the  lowest-order  func¬ 
tions  are  the  classical  aberrations  of  tilt,  focus,  coma,  etc.,  most  of  which 
are  known  to  have  large  values  in  phase  distortions  produced  by  turbulence 
and  thermal  blooming.  We  have  given  the  name  ’'ZEP-COAT"  to  a  multi¬ 
dither  modal  COAT  system  that  uses  Zernike  modes. 

Using  computer  simulation  codes,  we  have  determined  that  only  the 
7  polynomials  indicated  with  an  asterisk  in  Table  1  are  necessary  for  good 
correction  in  moderate  thermal  blooming.  In  determining  that  these 
7  polynomials  were  the  significant  ones,  a  moderate  blooming  case  was 
chosen  (^^^ansmit  ^  ^critical^*  desired  cw  correction  phasefront  was 

SL 

calculated  using  a  function -maximization  routine.  '  This  phasefront  was 

P.  Brown,  Jr.,  "Computer  Simulation  of  Adaptive  Optical  Systems,  ” 
Final  Report  on  Contract  N60921-74-C-0249,  September,  1975. 

J.  E.  Pearson,  et  al.  ,  ’’High  Power,  Closed  Loop,  Adaptive  System 
(HICLAS)  Study,"  Final  Report  on  Contract  N60921 -76-C-0608. 
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Table  1.  Zernike  Polynomials  Used  in  ZEP-COAT  Simulation 


n 

Z 

n 

Aberration  Name 

_ _ _ — - - — - 

Pn"'"' 

X 

Horizontal  tilt 

2 

2^ 

Y 

Vertical  tilt 

2 

*3 

2R^  -  1 

Refocus 

1.5 

>Jc4 

Astigmatism 

3 

5^ 

XY 

Astigmatism 

3 

X(3R^  -  2) 

Coma  (Y-symmetry) 

4 

7a 

Y(3R^  -  2) 

Coma  (X- symmetry) 

4 

8 

X(X^  -  3Y^) 

120  degree  (5th  order, 
Y-symmetry) 

4 

9^ 

Y(Y^  -  3X^) 

120  degree  (5th  order, 

X- symmetry) 

4 

*10 

Or'^  -  6R^  +  1 

Spherical  aberration 

2.  5 

11 

(X^  -  Y^)(4R^  -  3) 

Astigmatism  (5th  order) 

5 

12 

XY  (4R^  -  3) 

Astigmatism  (5th  order) 

5 

•A, 

"r 

13 

X(10R^  -  12R^  +  3) 

5th  order  coma  (Y-symmetry) 

6 

14^ 

YdOR**  -  12R^  +  3) 

5th  order  coma  (X-symmetry) 

6 

*15 

20R^  -  30R^  +  12R^  -  1 

5th-order  spherical  aberration 

3.  5 

where 

2  2 

X  =  -4=  =  R  cos  e,  Y  =  =  R  sin  6,  and  R  =  x  +  y  Si 

/2  /2 

^These  five  polynomials  are  not  u 

sed  for  compensating  only  thermal 

blooming  because  of  the  symmetries  in  the  blooming  distortions;  the 

transverse  wind  is  in  the  X-direction  (see  Ref*  11)*  They  are  used  for 
turbulence  compensation* 

'“'Principal  polynomials  needed  for  blooming  compensation* 

Relative  dither  amplitude  for  constant  effective  modulation  index. 

T1994 
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then  fitted  by  means  of  an  expansion  using  the  first  44  Zernike  polynomials. 
Table  2  lists  the  results  of  this  calculation  and  Figure  2  shows  the  uncor¬ 
rected  and  corrected  beam  profiles  calculated  for  the  conditions  shown  in 
Table  3.  Although  there  are  some  contributions  from  higher  order  poly¬ 
nomials,  their  amplitudes  are  either  less  than  1/10  of  the  smallest  ampli¬ 
tude  among  the  7  polynomials  noted  in  Tables  1  and  2,  or  are  of  too  high  an 
order  to  be  produced  by  a  37-element  mirror.  For  the  planned  experi¬ 
ments,  however,  we  will  implement  all  14  polynomials.  The  ZEP-COAT 
system  can  thus  be  used  for  turbulence  compensation  studies  if  so  desired. 

Implementing  a  ZEP-COAT  capability  into  the  DARPA/RADC  hard¬ 
ware  is  a  straightforward  procedure,  as  indicated  in  Figure  3.  A  ’’poly¬ 
nomial  generator”  is  merely  inserted  between  the  control  channels 
(synchronous  detectors  plus  low-pass  filters)  and  the  drivers  for  the 
deformable  mirror.  The  polynomial  generator  is  composed  of  several 
groups  of  weighting  matrices  that  distribute  the  error  signal  across  the 
mirror  in  the  correct  manner  to  produce  the  desired  polynomial  spatial 
mode.  Since  the  dither  and  correction  signals  will  be  put  onto  a  single 
mirror  for  our  experiments,  the  low-voltage  dither  and  correction  signals 
will  be  combined,  put  through  the  polynomial  box,  separated  and  separately 
amplified,  and  then  recombined  onto  the  mirror.  The  process  is  illustrated 
in  Figure  4. 

A  14-channel  polynomial  generator  is  nearly  complete;  37  mirror 
actuator  drivers,  each  capable  of  ±150  V  output,  and  a  dither/correction 
combiner  (dither  injector)  have  already  been  completed.  The  interconnec¬ 
tion  of  these  units  for  either  a  zonal  or  a  modal  control  system  is  shown  in 
Figure  5.  Because  of  the  excessive  delays  in  the  delivery  of  the  piezo¬ 
electric  transducer  (PZT)  material  for  the  deformable  mirror  (Section  B 
of  this  chapter),  the  technical  portion  of  the  contract  will  be  extended  in 
time,  but  at  no  increase  in  cost,  until  the  middle  of  November  1976  in 
order  to  complete  the  necessary  thermal  blooming  studies. 
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Figure  2,  Thermal  blooming  computer  simulation  data  for 
conditions  listed  in  Table  3.  (a)  Uncorrected 

focal-plane  beam  profile:  Strehl  ratio  =  0.  11. 

(b)  Corrected  beam  using  function-maximization 
routine;  Strehl  ratio  =  0,  34.  The  phase  correc¬ 
tion  used  here  is  used  to  generate  the  Zernike 
polynomial  amplitudes  in  Table  2. 
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Figure  4.  Zernike  polynomial  generator  schematic  for  the  DARPA/RADC  COAT  system.  The 
dither  and  correction  signals  are  combined,  put  through  the  polynomial  generator, 
separated,  amplified,  and  recombined. 
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Table  2.  Expansion  of  Thermal  Blooming  Phase  Correction  in  Terms  of 
Zernike  Polynomials  R^(p)  •  sin  m6  on  R^(p)  •  cos  m6 


Number 

m 

» 

Angular 

Function 

Radial  Polynomialj 

Aberration  Name 

Amplitude  in 
Blooming  Correction 
(radians) 

=:=! 

0 

2  ■ 

1 

2p^  -  I 

Refocus 

-1.  52 

*2 

0 

4 

1 

6p^  -  6p^  +  1 

Spherical  aberration 

-0.43 

«=3 

0 

6 

1 

20p^  -  SOp"^  r  12p^  -  I 

Spherical  aberration  (5th  order) 

-0.47 

4 

0 

8 

1 

70p®  -  140p^  +  OOp'*  -  20p^  r  1 

-0.  59 

1 

1 

sin 

P 

y-tilt 

0 

*6 

1 

1 

cos 

P 

x-tilt 

-10. 92 

*7 

1 

3 

sin 

3p^  -  2p 

coma 

0 

>:=8 

1 

3 

COS 

3p^  -  2p 

coma 

0.  19 

1 

5 

sin 

lOp^  -  12p^  +  3p 

coma  (5th  order) 

0 

*10 

1 

5 

cos 

lOp’  -  12p^  +  3p 

coma  (oth  order) 

0.36 

1 1 

1 

7 

sin 

35p"^  -  60p^  +  30p^  -  4p 

0 

12 

1 

7 

cos 

35p^  -  60p^  +  30p^  -  4p 

0,  64 

#13 

2 

2 

sin 

2 

P 

Astigmatism 

0 

#14 

2 

2 

COS 

2 

P 

Astigmatism 

3.95 

15 

2 

4 

sin 

4p'‘  -  3p^ 

Astigmatism  (5th  order) 

0 

*16 

2 

4 

cos 

4p^  -  3p2 

Astigmatism  (5th  order) 

0.  53 

17 

2 

6 

sin 

1  5p^  -  20p^  +  6p^ 

0 

18 

2 

6 

COS 

1  5p^  -  20p^  +  6p^ 

0.  26 

19 

2 

8 

sin 

S6p^  -  105p^  +  bOp"^  -  lOp^ 

0 

20 

2 

8 

cos 

56p^  -  105p^  +  60p^  -  lOp^ 

-0.  12 

*21 

3 

3 

sin 

3 

P 

120° 

0 

*22 

3 

3 

cos 

3 

P 

120° 

-0.  38 

23 

3 

5 

sin 

-  5  .3 

op  -  4p 

0 

24 

3 

5 

cos 

c  5  ,3 

5p  -  4p 

-0.  53 

25 

3 

7 

sin 

21p^  -  30p‘’  +  lOp^ 

0 

26 

3 

7 

cos 

21p^  -  SOp"*  +  lOp^ 

.0.62 

27 

4 

4 

sin 

4 

P 

0 

28 

4 

4 

cos 

4 

P 

0.27 

29 

4 

1  ^ 

sin 

6p^  -  5p^ 

0 

30 

4 

6 

cos 

6p^  -  5p^ 

‘0.71 

31 

4 

8 

sin 

28p®  -  42p^  +  15p^ 

0 

32 

4 

8 

cos 

28p®  -  42p^  t 

1.09 

33 

5 

5 

sin 

5 

P 

0 

34 

5 

5 

cos 

5 

P 

-0,56 

35 

5 

7 

sin 

Ip'^  -  6p’ 

0 

36 

5 

7 

cos 

■7  7  ,5 

7p  -  6p 

-0.  67 

37 

6 

6 

sin 

6 

P 

0 

38 

6 

6 

cos 

6 

P 

0.  58 

39 

6 

8 

sin 

8p®  -  7p^ 

0 

40 

6 

8 

cos 

8p^  .  7p^ 

0.  30 

41 

7 

7 

sin 

7 

P 

0 

42 

7 

7 

cos 

7 

P 

-0.  41 

43 

8 

8 

sin 

8 

P 

0 

44 

8 

8 

cos 

8 

P 

0.72 

*To  be  implemented  in  experimental  hardware. 


Table  3.  Assumed  Conditions  for  Thermal  Blooming  Calculations 


Range:  2  km 

Absorption  coefficient:  4  x 

ir.-4  -1 

10  m 

Wavelength:  10.6  pm 

Slew  rate:  10  mrad/sec 

Transmitter  diameter:  1  m, 
uniform  beam 

Transverse  wind  velocity: 

2  m/ sec  in  +X 
direction 

Resultant  uncorrected  Strehl  ratio:  0,  39 

T1996 


The  computer  simulations  being  used  to  analyze  ZEP-COAT 
performance  were  discussed  in  the  previous  contract  report.^  We  also 
reported  that  we  were  observing  an  instability  in  the  control  system.  We 
have  determined  that  the  instability  was  caused  by  an  improper  setting  of 
the  gain  and  modulation  index  for  the  channels.  We  have  now  optimized  both 
these  quantities  and  the  ZEP-COAT  servo  system  exhibits  stable  behavior 
under  all  conditions.  If  the  dither  signals  applied  to  the  mirror  are  repre¬ 
sented  by 


sinw^t)  Z^(r,  e)  ,  (3) 

n 

where  Zernike  polynomials,  the  the  relative  values  of 

for  a  5%  effective  modulation  index  are  given  in  the  last  column  of  Table  1. 
The  values  of  depend  on  the  order  of  the  polynomial  according  to 


n  +  1  ,  m  ^  0 

where  m  and  n  are  given  in  Table  2. 
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The  computer  simulations  have  also  been  used  to  investigate 
multiple  convergence  states  for  ZEP-COAT  systems,  that  is,  the  modal 
equivalents  of  2Nit  ambiguities  in  zonal  systems. 

The  2Ntt  ambiguity  problem  is  normally  associated  only  with  zonal 
multidither  systems.  An  analogous  problem  can  occur  with  a  modal  ZEP- 
COAT  system,  however.  That  is,  there  exist  local  maxima  states  where 
the  dither  control  system  can  stop  (obtain  a  zero  error  signal),  but  which 
are  nonoptimum  for  maximum  focal-plane  irradiance.  Taken  individually, 
the  two  most  obvious  cases  are  tilt  (the  servo  can  place  a  sidelobe  of  the 
Airy  pattern  on  a  glint)  and  focus. 

The  behavior  of  the  focus  control  can  be  understood  by  referring 
to  Figure  6.  Figure  6  illustrates  the  well-known  behavior  of  the  on-axis 
intensity  of  a  focused  beam:  as  the  focal  point  is  approached  from  either 
direction,  the  on-axis  intensity  goes  through  several  maxima  and  minima. 
Each  maxima  away  from  the  focal  plane  corresponds  to  a  different  Fresnel 
zone,  or  a  2N-it  phase  difference  between  rays  originating  from  the  center 
of  the  focusing  lens  and  the  outer  edge.  If  the  initial  defocus  is  at  a  point 
such  as  A  in  Figure  6,  the  ZEP-COAT  focus  control  will  drive  the  focus 
toward  the  maximum  on-axis  intensity.  If  the  starting  point  is  at  B,  how¬ 
ever,  the  multidither  servo  moves  toward  the  nearest  maximum  at  C  and 
a  nonoptimum  target  intensity  results. 

The  problem  can  also  exist  for  other  Zernike  functions.  Table  4 
lists  the  first  13  polynomials  along  with  the  location  of  the  first  secondary 
maximum  and  the  on-axis  irradiance  at  that  maximum.  At  this  point,  we 
have  only  pointed  out  the  potential  for  a  servo -ambiguity  for  ZEP-COAT. 
Whether  or  not  such  ambiguities  will  occur  in  actual  systems  when  they 
attempt  to  correct  for  arbitrary  or  random  errors  is  a  subject  for  future 
study. 

We  are  currently  setting  up  a  computer  code  which  will  determine 
the  relative  actuator  drives  necessary  to  produce  a  given  Zernike  poly¬ 
nomial  with  a  deformable  mirror.  This  code  will  be  used  to  set  the  values 
of  the  various  resistors  in  the  weighting  matrices  of  the  Zernike  polynomial 
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Figure  6.  Behavior  of  on-axis  intensity  of  a  focused  beam  showing 
secondary  maxima. 
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Table  4.  2Nir  Problem  with  Zernike  Polynomials 


Zernike  Polynomial  (Aberration) 

Approximate 
Location  of 
Secondary 
Maximum, 
(radians) 

On-Axis  Intensity 
at  Secondary 
Maximum  Relative 
to  that  of 

Unaberrated  Focus, 

1.  =  Pj(2p^  -  1)  [defocus] 

4. 5 

0.05  I 

0 

2.  4)  =  p2(6p'^  - 

5. 75 

0.  12  I 

0 

[spherical  aberration] 

3.  4)  =  -  3p  +  12p  -  1) 

8.  7 

0.  11  I 

0 

4.  4>  =  P4X  5.  4>  =  PgY 

5.  25 

0.017  I 

0 

[tilt] 

6.  4>  =  P6X(3p^  -  2) 

Tiny  one  at 

7.  5 

0.00012  I 

0 

[coma] 

7.  4>  =  P6Y(3p^  -  2) 

Next  at  12 

0.022  I 

0 

8.  4)  =  PgXllOp'^  -  12p^  +  3) 

9.  4>  =  PgYllOp'^  -  12p^  +  3) 

14 

0.012  I 

0 

10.  4>  = 

8.  5 

0.022  I 

0 

[astigmatism] 

11.  4,  =  Pii(X^  -  Y^)(4p^  -  3) 

13.  5 

0.03  I 

0 

12.  4>  =  Pi  2^^^^  “ 

None  for 
p.  <  25  rad 

— 

13.  4>  =  Pi3Y(Y^  -  3X^) 

1 

T1997 
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generator.  The  best  fit  to  a  particular  is  a  function  of  the  mirror 
influence  function,  the  number  of  actuators,  and  the  interactuator  spacing. 
Proper  representation  of  a  is  necessary  in  order  to  maintain  the 
orthogonality  of  the  COAT-controlled  functions.  If  the  mirror  cannot  accu¬ 
rately  produce  the  Zernike  functions,  the  functions  that  are  produced  will 
not  be  orthogonal  and  the  COAT  control  channels  will  be  coupled  together. 
The  amount  of  the  coupling  will  be  a  function  of  the  rms  difference  between 
the  mirror-produced  function  and  the  desired  Z^, 

Another  potential  source  of  coupling  with  a  modal  system  such  as 
ZEP-COAT  occurs  through  the  amplitude  distribution  across  the  laser 
transmitting  aperture.  The  polynomials  listed  in  Table  1  are  orthogonal 
over  a  uniform  unit  circle.  For  a  nonuniform  intensity  distribution,  they 
will  not  be  orthogonal.  If  the  intensity  distribution  is  known,  a  new  poly¬ 
nomial  set,  related  to  the  Zernike  set,  can  be  found  that  is  orthogonal.  If 
the  polynomials  are  not  orthogonal,  the  result  is  again  a  coupling  of  the 
COAT  servo  channels. 

Just  how  much  channel-channel  coupling  a  ZEP-COAT  system  can 

2  3 

tolerate  is  not  known.  Previous  studies  of  zonal  systems  *  have  shown 
that  20  to  30  percent  coupling  between  channels  can  be  tolerated  in  a  sine/ 
cosine  dither  system  without  noticeable  degradation  of  performance.  To 
our  knowledge,  this  important  problem  of  channel-channel  coupling  in  a 
modal  control  system  has  not  been  satisfactorily  investigated. 

Deformable  Mirror 

For  the  experimental  studies  with  ZEP-COAT,  a  37-actuator 
deformable  mirror  is  necessary,  and  one  is  being  built  to  replace  the 
beam-splitter/phase-shifter  array  (called  the  ’’phasor  matrix*’)  that  has 
been  successfully  used  up  to  now  with  the  DARPA/RADC  COAT  system. 

The  mirror  will  perform  the  two  functions  of  phase  dithering  (tagging)  and 
phase  correction  required  in  any  type  of  multidither  COAT  system.  Per¬ 
formance  of  both  these  functions  with  a  single  element  places  heavy  and 
conflicting  demands  on  the  devices.  The  phase  correction  function  requires 
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a  large-amplitude  phase  excursion  (about  ±2  optical  wavelengths)  at 
relatively  low  frequencies  (up  to  about  1  kHz),  On  the  other  hand,  the 
dither  function  requires  low-amplitude  excursions  (±k/10)  at  much  higher 
frequencies  (10  to  30  kHz).  Xo  date,  no  continuous  — surface,  deformable 
mirror  device  has  been  constructed  that  can  accomplish  both  of  these  func¬ 
tions  simultaneously,  particularly  with  the  low  drive  voltages  (±125  V) 
available  from  the  DARPA/RADC  COAT  system. 

The  mirror  being  built  on  this  program  is  the  first  of  a  kind,  and 
as  such  will  represent  the  state  of  the  art  for  uncooled  deformable  mirrors 
of  its  type.  Because  such  a  mirror  had  not  previously  been  built,  or  even 
designed,  an  extensive  design  effort  was  initiated.  This  design  has  been 
funded  only  in  part  by  this  contract,  the  remainder  of  the  work  being  per¬ 
formed  as  part  of  the  Hughes  Aircraft  IR&D  program. 

The  design  and  analysis  effort  started  with  Hughes  proprietary 
designs  for  multielement  deformable  mirrors  and  for  external- spring,  self- 
contained,  piezoelectrically-driven  actuator  cells.  The  basic  design  and 
performance  goals  are  discussed  in  an  earlier  report. 

A  schematic  of  the  mirror  is  shown  in  Figure  7.  The  entire  mirror 
body  is  made  of  beryllium.  The  37  actuators  are  placed  on  a  circular  array 
with  a  minimum  actuator  spacing  of  0.  5  50  in.  The  active  mirror  surface 
has  a  diameter  of  3.  980  in.  The  faceplate  is  initially  0.  150-in.  -thick 
beryllium  optical  stock,  and  is  butt-brazed  to  each  of  the  37  actuators  and 
to  the  rim  of  the  primary  backup.  The  faceplate  will  be  polished  and  coated 
with  a  protective  silver  overcoat.  The  actuator  cells  are  brazed  to  the  pri¬ 
mary  backup  at  the  same  time  as  the  faceplate  braze. 

Each  piezoelectric  (PZT)  transducer  stack  is  electrically  grounded 
at  the  mirror  faceplate  and  the  driver  voltage  is  applied  through  the  elec¬ 
trode  at  the  other  end  of  the  PZT  stack.  With  this  arrangement,  the  trans¬ 
ducer  stack  is  wired  internally,  thus  eliminating  the  complication  of  feedback 
wires  out  from  the  stack.  The  only  wire  to  each  actuator  is  soldered  con¬ 
veniently  to  the  electrode  in  a  manner  which  directs  the  wire  straight  out  of 
the  cell,  with  no  bends  that  require  additional  volume.  This  electrical 
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arrangement  is  compatible  with  the  electronics  of  the  DARPA/RADC 
system,  which  has  a  common  ground  for  all  channel  outputs.  All  of  the 
electrical  outputs  are  brought  to  a  41 -pin  electrical  connector  that  is 
attached  to  the  dust  cover  at  the  back  of  the  mirror  (see  Figure  7). 

The  mirror  has  been  fabricated  and  is  currently  being  polished, 
although  only  18  of  the  required  37  PZT  stacks  are  presently  available.  A 
photograph  of  the  unit  is  shown  in  Figure  8.  An  analysis  performed  on  this 
mirror  predicts  that  the  lowest  structural  resonance  will  occur  at  12  to 
15  kHz,  and  appears  to  be  a  shear  mode  of  the  back-plate  structure.  The 
measured  response  characteristics  for  the  central  actuator  (before  polish¬ 
ing)  is  shown  in  Figure  9,  with  all  other  actuators  undriven,  but  with  pre-  ^ 

loads  on  them  (dummy  blocks  were  loaded  into  the  outer  ring).  The 
response  curves  were  obtained  by  placing  the  mirror  in  a  Twyman-Green 
interferometer  and  looking  at  the  surface  motion  of  the  central  actuator. 

The  response  of  the  other  actuators  is  similar  except  for  some  of  the 
details  of  the  fluctuations  beyond  20  kHz.  The  phase  of  the  surface  response 
is  measured  relative  to  the  sinusoidal  drive  voltage. 

The  mirror  response  is  reasonably  flat  up  to  about  20  kHz,  where  a 
rather  complicated  resonance  behavior  beings.  Using  a  holographic 
interferometer,  we  conclude  that  the  fluctuations  above  20  kHz  are  a  result 
of  excitation  of  complicated  actuator/faceplate  vibrational  modes.  The 
first  two  major  resonances  of  this  type  were  calculated  to  occur  at  23  kHz 
and  40.  5  kHz,  in  good  agreement  with  the  first  two  major  peaks  in  the 
amplitude  plot  in  Figure  9. 

Note  that  the  first  resonance  activity  occurs  in  the  neighborhood  of 
15  kHz,  which  is  in  good  agreement  with  the  structural  analysis.  The 
relatively  small  phase  fluctuations  indicate  that  the  mirror  is  useable  out 
to  at  least  30  kHz.  Even  a  higher  frequency  response  may  be  possible  by 
filling  the  interior  of  the  unit  with  oil,  and  by  using  lower-amplitude  drive 
voltages  at  the  higher  frequencies  (as  will  be  done  for  dithers)  to  reduce 
the  drive  to  backplate  modes.  At  this  time  we  do  not  fully  understand  the 
reasons  for  the  enhanced  amplitude  response  above  20  kHz. 
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Figure  8. 


Photograph  of  37-element  DARPA/RADC  deform¬ 
able  mirror  prior  to  final  polishing. 
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In  addition  to  the  frequency  response  measurements  shown  in 
Figure  9,  we  have  measured  the  influence  function  of  an  actuator  in  each 
of  the  3  rings  plus  the  central  actuator.  Profiles  through  the  center  of 
these  4  actuators  are  shown  in  Figure  lOo  Each  actuator  can  produce 
0.28  \xm  of  surface  motion  for  150  V  of  drive.  As  can  be  seen,  the  influ¬ 
ence  function  depends  on  the  distance  of  the  actuator  from  the  clamped 
edge  of  the  mirror;  the  closer  to  the  edge,  the  more  skewed  the  surface 
profile. 

The  mirror  surface  is  measured  using  a  simple,  but  highly  sensi¬ 
tive  and  accurate  technique.  A  sinusoidal  drive  at  500  to  1000  Hz  is  first 
applied  to  the  appropriate  actuator  (or  group  of  actuators).  About  30  V  rms 
is  used,  which  produces  a  peak  surface  deflection  of  about  0.079  iJim.  A 
phonograph  needle  is  then  drawn  across  the  surface  and  the  resulting  out¬ 
put  is  observed  on  an  rms  voltmeter  as  a  function  of  the  needle  location. 

The  observed  profiles  are  repeatable  to  within  a  few  percent.  The  tech¬ 
nique  is  particularly  useful  because  it  gives  the  actuator  influence  function 
independent  of  the  undriven  mirror  surface  shape.  For  use  on  polished 
mirrors,  the  needle  tip  is  dipped  in  epoxy  to  form  a  nonabrasive  rounded 
surface.  The  epoxy  produces  very  little  loss  in  sensitivity. 

As  indicated  in  Figure  10(a),  the  central  actuator  influence  function 
is  closely  approximated  by  the  "subgaussian”  function, 

IN(r)  =  exp  (-3.  622r^  *  ^)  ,  (5) 

where  r  is  the  radial  distance  in  inches  from  the  actuator  center.  Fig¬ 
ure  11  shows  a  3-dimensional  plot  of  this  function  and  Figure  12  is  a 
contour-line  plot  of  this  profile.  Notice  that  the  profile  is  reasonably 
symmetric,  but  is  slightly  higher  at  the  location  of  the  actuator  marked 
(a)  than  at  other  actuator  locations. 

If  we  define  a  mechanical  coupling,  C^,  as 

=  IN(S)  ,  (6) 


31 


NORMALIZED  SURFACE  NORMALIZED  SURFACE 

DEFLECTION  DISPLACEMENT 


CENTER  OF 
NEAREST 
ACTUATOR 


CENTER  OF 
NEAREST 
ACTUATOR 


exp(-3.622rl-5)- 


-1.0  -0.8  -0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.8  1.0 

r=  RADIAL  DISTANCE,  in. 


LU 

I  TOWARD  EDGE 

o 

< 

OFMIRROR  /"\ 

u. 

QCZ 

0.8 

DO 

/  \ 

QO 

0.6 

/  \ 

U.1UJ 

tsj-J 

1  /  \  1 

i- 

SO 

0.4 

s 

a: 

0.2 

O 

z 

0 

]  1  1  1  1  > _ 1 _ 1^  ^ 

-0.8  -0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.8  1.0 

r=  RADIAL  DISTANCE,  in. 

1.0  h  TOWARD  EDGE 
OF  MIRROR 


CLEAR 
APERTURE 
EDGE  i 


qH  ocL  OUTER 
So  0.6  h  edge  OF 


-0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.8  1 .0 

r  =  RADIAL  DISTANCE,  in. 


-0.8  -0.^0.4  -0,2  0  0.2  0.4  0.6  0.8  1 .0 

r=  RADIAL  DISTANCE,  in. 


Figure  10.  Influence  function  profiles  of  beryllium  mirror. 

(a)  Center  actuator,  (b)  First  ring  actuator. 

(c)  Second  ring  actuator,  (d)  Third  (outermost)  ring 
actuator.  An  empirical  curve  fit  to  the  central  actu¬ 
ator  profile  is  shown  as  a  dotted  line  in  (a);  the  curve 
is  exp  (-3.  122  r^‘  ^),  where  r  is  the  radial  distance 
from  the  profile  peak  (actuator  center)  in  inches.  Tl 
inter-actuator  center -center  spacing  is  0.  55  in. 


Figure  11. 


Three-dimensional  view  of  centeral-actuator  influ¬ 
ence  function  for  RADC  beryllium  mirror. 
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where  S  is  the  interactuator  center-center  spacing,  then  Cj^  =  0.  23  for  this 
RADC  beryllium  mirror. 

The  more  important  quantity,  however,  is  the  servo  coupling,  • 

O'Meara^  has  defined  this  coupling  as  the  ratio  of  the  error  signal  induced 
in  a  given  channel  by  the  displacement  of  a  single  neighboring  actuator  to 
that  induced  by  a  comparable  displacement  of  the  actuator  associated  with 
the  given  channel.  When  the  influence  function  is  of  the  form 


IN(r)  =  exp 


(7) 


2  2  2 

where  r  =  x  +  y  is  the  squared  distance  from  the  actuator  center,  then 
the  servo  coupling  is"^ 


C 

sn 


2Trr{2/n)  n’ 


where  n  is  any  number,  r(x)  is  the  gamma  function,  and 


I 

n 


(X  - 


P)^  + 


in/2 


dxdy. 


(9) 


The  quantity  p  is  given  by 


(-inC  ) 
m 


1/n 


(10) 


Reference  4  gives  the  values  of  C  for  various  choices  of  C  and  n.  For 

sn  m 

the  RADC  mirror,  n  =  1.  5  and  C  =0.  23,  giving  C  =0.  70.  For  C  =0.  20, 

m  ’  ®  sn  m  ’ 

C  =  0.67  and  if  C  =  0.  15,  C  =  0.58.  All  of  these  values  are  very 
sn  m  ’  sn  ’ 

substantial  couplings,  which  may  severely  limit  the  usable  COAT  servo 
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bandwidth  when  this  mirror  is  employed.  This  is  a  point  which  we  will 
have  to  verify  experimentally. 

The  only  way  to  reduce  the  servo  coupling  once  n  is  fixed,  Eq.  (3), 

is  to  narrow  the  influence  function,  i.  e. ,  reduce  Studies  conducted 

on  the  Hughes  IR&D  program,  however,  show  that  n  =  2.  0  minimizes 

C  for  a  fixed  value  of  C  .  On  the  other  hand,  reducing  C  will  increase 
sn  m  ^ 

the  inter  actuator  ripple  (surface  ripple  when  all  actuators  have  equal 
excitation).  The  tradeoff  between  ripple  and  as  a  function  of  is 

shown  in  Figure  13.  The  ripple  data  are  obtained  by  calculating  the  surface 
produced  by  a  99-actuator,  hexagonal-array  mirror.  Each  actuator  is 
given  a  unit  amplitude  deflection  and  all  actuators  are  assumed  to  have  the 
same  influence  function.  The  resulting  surface  is  assumed  to  be  a  linear 
superposition  of  all  the  individual  actuator  influence  functions.  This  is  an 
approximation,  but  one  which  should  be  reasonably  good  as  long  as  the 
actuator  restoring  springs  dominate  the  interactuator  surface  ’’springs*’ 
caused  by  the  faceplate  stiffness  (this  is  expected  to  be  the  case  in  the  RADC 
mirror). 

Note  in  Figure  13  that  curves  (b)  and  (c)  for  the  ripple  cross  over. 

This  implies  that,  for  a  given  C^,  there  is  a  value  of  n  in  Eq,  (7)  that 

minimizes  the  ripple.  This  optimum  exponent  is  plotted  in  Figure  14  as 

a  function  of  C  .  For  large  C  ,  the  value  of  n  .  approaches  2.  0,  a  gaus- 
m  m  opi 

sian.  For  all  cases,  however,  ^  RADC  mirror  is  expected 

to  have  a  larger  surface  ripple  when  excited  than  a  mirror  with  a  gaussian 
influence  function. 

In  fact,  based  on  the  data  in  Figures  13  and  14,  it  appears  that  the 
influence  function  on  the  RADC  mirror  is  far  from  optimum.  A  value  of  n 
between  2  and  2.  5  is  desirable,  with  i  10%  so  that  the  servo  coupling 
will  be  less  than  30%. 

The  experimentally  observed  surface  profile  when  two  adjacent 
actuators  are  driven  with  the  same  voltage  is  shown  inFigure  15.  For 
reference,  the  individual  actuator  profiles  are  also  shown.  The  peak-to- 
peak  surface  ripple  is  4.  9%  of  the  surface  deflection  at  one  actuator  loca¬ 
tion.  Figure  13  predicts  8%  ripple  if  aU  actuators  are  energized.  A 
simple  superposition  of  2  actuators  with  the  same  influence  function. 
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Figure  13.  Servo  coupling  C  (eq.  (8))  and  peak-peak  ripple 
as  a  function  of  deformable  mirror  mechanical 
coupling,  Cni  (eq.  (6)).  Three  influence  function 
cases  are  shown,  a  "subgaus sian”  (n  =  1.5),  a 
gaussion  (n  =  2.  0),  and  a  ’’supergaussion''  (n  =  2.  5). 
Further  data  are  available  in  Ref.  4. 
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=  OPTIMUM  EXPONENT 


5554-2 


Figure  14.  Optimum  value  of  exponent,  n,  in  eq.  (3)  in  order  to 
produce  minimum  surface  ripple. 
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Figure  15,  Experimentally  observed  surface  profile  with  two 

adjacent  actuators  on  DARPA/RADC  mirror  driven 
with  equal  voltage  {^30  V  rms  at  500  Hz).  The 
individual  actuator  influence  functions  that  produce 
the  resultant  surface  are  also  shown. 
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IN  =  exp  [-  (r/r^)^*  ^],  and  23%  coupling,  predicts  a  3.  3%  ripple.  The 
difference  between  this  value  and  the  observed  4.  9%  can  be  attributed  to 
the  differences  between  the  two  influence  functions  in  Figure  14;  they  are 
not  identical. 

C.  Irradiance  Tailoring  Studies 

The  previous  contract  report^  presented  some  results  on  efforts  to 
find  a  fixed  transmitter  irradiance  distribution  that  will  minimize  thermal 
blooming  effects.  The  basic  conclusion  presented  there  was  that  special  laser 
modes  and  truncated  distributions  have  little  effect  on  the  strength  of  bloom¬ 
ing  at  a  fixed  transmitted  power.  In  order  to  minimize  blooming,  the  trans¬ 
mitter  aperture  should  be  made  as  large  as  possible  and  the  irradiance 
distribution  should  be  made  as  uniform  as  possible. 

Since  the  time  the  first  report  was  written,  we  have  found  an  error 
in  the  calculations  that  produced  Figures  7  and  10  of  that  report.  For¬ 
tunately,  the  basic  conclusion  remains  unchanged  even  though  the  quantita¬ 
tive  output  of  the  calculations  is  different.  The  corrected  work  has  been 
accepted  for  publication  in  the  Journal  of  the  Optical  Society  of  America, 
to  appear  in  early  1977.  Since  the  corrected  results  appear  in  that  article, 
we  present  it  here  in  its  entirety  as  an  Appendix. 
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III.  PLANS  FOR  THE  NEXT  CONTRACT  PERIOD 

The  next  phase  of  the  contract  will  be  devoted  almost  exclusively  to 
characterizing  the  RADC  beryllium  mirror.  The  faceplate  will  be  thinned 
to  0.  125  in.  to  reduce  the  mechanical  coupling,  and  the  unit  will  be  polished 
with  19  actuators  in  place.  All  of  the  data  required  by  Amendment  No,  1 
to  the  contract  will  be  obtained.  These  data  will  form  the  substance  of  the 
3rd  interim  report  of  the  contract. 

In  addition  to  completing  the  construction  and  characterization  of  the 

mirror,  we  will  also  complete  the  Zernike  polynomial  generator  and  set  up 

the  thermal  blooming  optical  system.  Our  present  plans  are  to  use  a 

flowing -liquid  blooming  medium  for  two  reasons:  (1)  difficulties  and  safety 

hazards  with  the  flowing-gas  cell;  (2)  ease  of  producing  substantial  blooming 

6)  —  8 

with  the  limited  optical  power  available.  Previous  work  has  indicated  no 
signficant  differences  between  thermal  blooming  produced  in  a  gaseous  or  in 
a  liquid  medium. 
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APPENDIX 


Propagation  of  Laser  Beams  Having  an  On- Axis  Null 
in  the  Presence  of  Thermal  Blooming 
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PROPAGATION  OF  LASER  BEAMS  HAVING  AN  ON- AXIS  NULL 
IN  THE  PRESENCE  OF  THERMAL  BLOOMING* 

J.  E.  Pearson,  C.  Yeht,  and  W.  P.  Brown,  Jr. 

Hughes  Research  Laboratories 
3011  Malibu  Canyon  Road 
Malibu,  California  90265 

ABSTRACT 

The  propagation  of  focused  beams  having  an  on-axis  irradiance 
null  is  considered  in  the  presence  of  thermal  blooming*  Two  cases  are 
treated:  (a)  beam  profiles  that  have  an  irradiance  zero  in  the  beam  center 
at  the  focal  plane  as  well  as  the  transmitter  aperture;  (b)  beam  profiles  that 
have  an  on-axis  irradiance  null  only  at  the  transmitter.  It  is  demonstrated 
that  none  of  the  beam  profiles  considered  in  case  (a)  has  a  meaningful 
advantage  over  a  Gaussian  beam  profile.  Some  of  the  case  (b)  profiles  do 
produce  a  larger  bloomed  irradiance  in  the  focal  plane,  particularly  when 
the  initial  intensity  distribution  is  very  uniform  in  its  non- zero  regions. 
Addition  of  a  simple  central  obscuration  to  a  ’’filled’*  irradiance  distribution 
is  found  to  have  no  advantage,  however,  for  the  cases  considered. 


*To  be  published  in  the  Journal  of  the  Optical  Society  of  America.  This 
work  sponsored  in  part  by  the  Air  Force  Systems  Command’s  Rome  Air 
Development  Center,  Griffiss  AFB,  NY. 

t  Permanent  address:  Electrical  Engineering  Department, 

University  of  California  at  Los  Angeles 
Los  Angeles,  California  90024 
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I,  INTRODUCTION 

When  an  optical  laser  beam  is  propagated  through  an  absorbing 
medium  with  a  negative  refractive  index  temperature  coefficient  (dn/dT  <  0), 
the  phenomenon  of  thermal  blooming  ’  occurs.  In  the  presence  of  thermal 
blooming,  the  propagation  medium  acts  like  a  negative  lens  that  defocuses, 
spreads,  and  distorts  the  optical  beam.  This  power -dependent  phenomenon 
can  be  particularly  harmful  when  the  optical  system  is  attempting  to  produce 
maximum  focal  plane  irradiance. 

It  has  been  suggested  that  the  use  of  a  transmitted  beam  that  has  an 
annular  intensity  profile  may  provide  higher  focal  plane  irradiance  in  the 
presence  of  thermal  blooming^  than  that  produced  by  a  Gaussian  beam.  The 
physical  reasoning  behind  this  suggestion  is  straightforward.  A  Gaussian 
intensity  profile  produces  a  quadratic  refractive  index  variation  in  the 
medium  that  is  minimum  at  the  beam  center  where  most  of  the  optical  power 
occurs.  With  an  annular  beam  that  has  an  on-axis  null,  the  negative  lens 
induced  in  the  medium  has  a  larger  refractive  index  on-axis  than  off-axis. 
The  result  is  a  positive  lens  for  light  that  is  close  to  the  beam  axis,  and 
this  lens  tends  to  counteract  some  of  the  off-axis  negative  lens  beam  spread¬ 
ing  as  the  optical  beam  propagates. 

This  article  presents  the  results  of  an  investigation  aimed  at 
quantifying  how  much  thermal  blooming  distortions  are  reduced  by  the  use 
of  laser  beam  that  have  an  on-axis  intensity  zero.  Such  beams  can  be 
produced  either  by  apodization  (truncation),  or  by  use  of  laser  resonators 
that  have  annular  gain  profiles,^  or  by  many  unstable  resonator  con¬ 
figurations  such  as  those  used  in  many  high  energy  laser  devices.  The 
investigations  were  performed  by  computer  simulation  of  the  propagation 
problem.  The  computer  code^  solves  the  coupled  differential  equations 
that  describe  both  the  medium  dynamics  and  the  diffraction-propagation 
of  a  focused  laser  beam.  Convection-dominated  heat  transfer  produced  by 
a  transverse  wind  is  assumed  and  the  quantity  of  interest  is  the  peak 
focal  plane  irradiance. 

Two  types  of  annular  or  "on- axis -null"  intensity  profiles  are  treated. 
Section  II  discusses  those  that  have  an  intensity  zero  in  the  beam  center. 
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both  in  the  focal  plane  (far -field)  as  well  as  at  the  transmitter  (near -field) 
and  thus  are  free- space  propagation  modeSo  Section  III  discusses  those 
profiles  that  have  a  zero  in  the  beam  center  only  at  the  initial  transmitter 
plane.  Due  to  diffraction,  the  initial  on-axis  null  fills  in  as  these  beams 
propagate  so  that  the  on-axis  irradiance  is  not  zero  at  the  focal  plane.  An 
appropriate  choice  of  the  initial  profile  can  indeed  increase  the  maximum 
focal  plane  irradiance.  The  standard  of  comparison  in  all  cases  is  the 
focal  plane  irradiance  produced  by  a  Gaussian  beam  that  is  truncated  at  the 
10%  intensity  radius. 

n.  "ON-AXIS-NULL”  IRRADIANCE  PROFILE  AT  THE  TRANSMITTING 
APERTURE  (NEAR-FIELD)  AND  AT  THE  FOCAL  PLANE 
(FAR -FIELD) 

In  this  section,  we  consider  cases  where  the  irradiance  null  in  the 

middle  of  the  beam  is  retained  as  the  focused  beam  is  propagated  to  the 

target.  To  first  demonstrate  that  an  appropriately  chosen  transmitter 

irradiance  distribution  can  indeed  be  focused  to  produce  an  on-axis-null 

focal-plane  irradiance  distribution,  we  start  with  a  transmitter  aperture 
4 

field  given  by 


E(r  Bj)  =  cos(mej),  (m  =  0,  1,  2  .  .  . . )  »  (1) 

where  A^  is  the  maximum  field  magnitude,  i^^ty  be  a  continuous  or 

a  discontinuous  function  of  r^,  and  (r^,  0^)  are  the  transverse  polar  coor¬ 
dinates  of  the  transmitting  aperture.  From  Fresnel  diffraction  formulae, 
we  obtain  the  field  u  at  the  focal  plane: 


-oo 
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where  z  is  the  focal  distance,  X.  the  free- space  wavelength,  k  =  2ir/\, 
are  the  transverse  rectangular  coordinates  at  the  focal  plane,  and  (Xj,yj) 
are  the  transverse  rectangular  coordinates  at  the  transmitting  aperture 
plane.  Changing  Eq.  (2)  into  polar  coordinates  and  carrying  out  the 
angular  integration,  we  find 
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The  on- axis  irradiance  occurs  at  r^  =  0  and  is  given  by 


u(0,  z) 


[Ztt  cos(m0^)T 

J 


(0)  . 

m' 


(5) 


Since  J  (0)  =  0  when  m  0,  Eq.  (4)  indicates  that  1=0  for  m  0.  An 
on-axis -null  focal  plane  irradiance  is  thus  obtained  for  all  m  ^  0.  When 
m  =  0,  however,  the  on- axis  focal  plane  irradiance  may  take  on  any  value, 
including  zero,  depending  upon  the  initial  radial  dependence  of  u^ 
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The  particular  transmitter  irradiance  profiles  that  we  have  chosen 
to  study  are  given  in  Eqs.  (6)  and  (7)  and  are  illustrated  in  Fig,  1.  The 
truncated  Gaussian  beam  we  are  using  for  comparison  purposes  is  also 
shown  in  Fig.  1  and  defined  in  Eq.  (8). 


2/  2 
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For  these  beams,  the  first  two  of  which  can  be  produced  by  coaxial- 

4 

resonator,  annular-gain  lasers,  is  the  e-folding  intensity  radius 

of  the  Gaussian  beam  in  Eq.  (8),  p  is  a  constant  less  than  or  equal  to  1, 
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Figure  1.  Three-dimensional  plots  of  transmitter  irradiance  profiles 
that  retain  their  shape  for  free- space  propagation.  The 
truncated  Gaussian  beam  is  used  as  a  comparison  reference  for 
all  other  beams  considered.  The  two  "on-axis-null"  profiles 
shown  here  also  have  an  on-axis  null  in  the  focal-plane. 


is  the  adjustable  maximum  intensity  radius,  and  and  a^  are  the 
inner  and  outer  radii  of  a  coaxial  laser  resonator.  The  beams  in  Eqs,  (6) 
and  (7)  are  assumed  oriented  in  such  a  way  that  the  =  tt/2  axis  is 
parallel  to  the  transverse  wind  direction.  This  orientation  produces 
minimum  blooming  since  the  heated  air  from  one  main  transmitter  lobe 
does  not  pass  through  the  region  heated  by  the  other  main  lobe.  The  field 
amplitudes  A^,  and  A^  are  chosen  so  that  the  total  transmitted  power 

is  constant.  The  expressions  used  to  obtain  these  amplitudes  are  derived 
in  the  appendix. 

Figure  2  compares  the  initial  transmitter  irradiance  distributions 
of  the  three  beams  to  the  focal-plane  irradiance  distributions  for  a  repre¬ 
sentative  transmitter  power.  The  parameters  necessary  to  specify  the 
blooming  strength  are  given  in  the  appendix.  For  all  the  computational 
results  presented  here,  we  have  chosen  values  for  the  various  constants  to 
correspond  to  the  obscuration  ratio  in  some  laser  pointer /tracker  systems 
now  in  use,  and  to  give  the  most  uniform  initial  intensity  distribution;  we  have 
found  that  thermal  blooming  is  minimized  for  uniform  initial  intensity 
distribution;  we  have  found  that  thermal  blooming  is  minimized  for  uniform 
initial  intensity  profiles.  The  values  used  in  our  calculations  are  ^ 

=  1.  52,  a  /b  =  7,  r  /a  =  2/3,  and  p  =  1,  Note  that  in  the  absence  of 
blooming  effects,  the  focal  plane  distributions  would  look  like  the  transmitter 
distributions  for  each  case. 

The  peak  focal- plane  irradiance  as  a  function  of  total  transmitted 
laser  power  for  these  three  beam  cases  is  given  in  Fig.  3.  Both  the 
irradiance  and  the  transmitter  power  are  plotted  in  normalized  units  as 
discussed  in  the  appendix.  The  maximum  focal -plane  irradiance  for 
Case  (Il-b)  is  a  factor  of  1.  2  larger  than  that  obtained  with  the  truncated 
Gaussian  beam.  The  higher  irradiance  is  obtained  by  increasing  the  total 
laser  power  at  the  aperture  plane  by  a  factor  of  2.  75.  Case  (Il-a)  yields 
lower  peak  focal-plane  irradiance  even  with  higher  total  laser  power. 

Hence,  no  significant  improvement  is  obtained  with  transmitting  intensity 
profiles  tailored  according  to  cases  (Il-a)  or  (Il-b). 
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Pseudo-gray  scale  plots  of  transmitter  and  focal-plane 
irradiance  distributions  for  the  three  beams  in  Fig.  1.  A 
moderate  amount  of  thermal  blooming  is  present  in  the  pro 
pagation  path. 
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III.  ”ON-AXIS-NULL"  INTENSITY  PROFILE  ONLY  AT  THE 
TRANSMITTING  APERTURE  (NEAR-FIELD) 

We  now  treat  cases  for  which  the  on-axis  intensity  null  occurs  only 
at  the  transmitting  aperture.  In  other  words,  under  thermal  blooming 
conditions,  the  cooler  region  in  the  middle  of  the  beam  does  not  exist 
throughout  the  focused  propagation  path  to  the  focal  plane.  Four  specific 
aperture  intensity  profiles  are  considered  and  compared  to  the  truncated 
Gaussian  in  Eq.  (8): 


Case 

(Ill-a) 


(9) 


Case 

(Ill-b) 


0 


,  for  b  <  r.  <  a 
2  o  1  o 

P 

,  otherwise 


(10) 


Case  |u  -  A? 
(III-c) 


=  0 


,  for  b  <  r .  <  a 
*  o  1  o 

,  otherwise 


(11) 


Case 

(Ill-d) 


2,  2 

I'Po 


=  0 


,  for  b  <  r,  <  a  ; 
*  o  1  o 

otherwise  . 
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The  beam  amplitudes  A^,  Ag,  A^  are  chosen  to  keep  the  total 
transmitter  power  constant.  The  calculations  for  these  amplitudes  are 
given  in  the  appendix.  Computer-generated  plots  of  the  initial  beam 
pforiles  are  illustrated  in  Fig.  4. 

The  transmitter  aperture  intensity  profiles  as  well  as  the  corre¬ 
sponding  focal  plane  intensity  profiles  for  these  cases  are  displayed  in 
Fig.  5  for  a  representative  laser  power  level.  Figure  6  compares  the 
normalized  peak  focal  plane  ir radiance  as  a  function  of  normalized  total 
input  transmitting  power  for  these  four  cases  to  that  for  the  truncated 
Gaussian  case.  Unlike  the  cases  treated  in  Section  II,  significantly  larger 
peak  focal  plane  irradiance  is  achieved  for  two  of  the  four  cases  considered 
as  compared  with  that  for  the  truncated  solid  Gaussian.  Case  (Ill-d)  yields 
a  peak  focal  plane  irradiance  which  is  the  same  as  that  for  the  truncated 
Gaussian  beam. 

The  intensity  profile  for  Case  (Ill-a)  is  obtained  by  combining  the 
TEMqq  and  TEM|'q  modes  of  a  cylindrical  cavity  oscillating  in  phase 
opposition  (Laguerre-Gaussian  modes^).  The  peak  focal  plane  irradiance 
achieved  with  this  profile  is  almost  twice  that  of  the  truncated  solid  Gaussian 
beam  with  a  total  laser  power  roughly  2.  5  times  that  for  the  truncated  solid 
Gaussian  beam.  Other  ways  to  increase  focal  plane  irradiance  by  this 

g 

amount  such  as  coherent  optical  adaptive  techniques  (COAT)  may  be  more 
practical  and  efficient  if  this  type  of  mode  is  difficult  to  produce  in  practice 
in  high  power  lasers. 

In  cases  (Ill-b)  and  (III-c),  the  peak  focal  plane  irradiance  is  also 
about  twice  the  value  for  the  truncated  solid  Gaussian  beam,  but  as  shown 
in  Fig.  6,  this  increase  is  obtained  at  about  90%  of  the  total  input  laser 
power  of  the  truncated  solid  Gaussian  beam.  A  very  uniform  intensity 
distribution  is  thus  desirable  for  reducing  blooming,  but  may  be  very 
difficult  to  obtain  in  practice. 

IV.  CONCLUSION 

These  results  lead  us  to  conclude  that  many  special  laser  modes  or 
intensity  distributions  that  have  an  on-axis  null  or  minimum  do  not  offer 
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CASE  m-b 


Figure  4. 


Three-dimensional  plots  of  transmitter  irradiance  profiles  that 
have  an  on-axis  null  initially  that  fills  in  as  the  beams  pro¬ 
pagate.  Case  Ill-a  is  a  mixture  of  cylindrical  modes  TEMqo 
and  TEMm  oscillating  in  phase  opposition.  Case  Ill-b  is 
similar  to  Case  ll-b,  except  |  u  |  ^  has  no  0 -dependence. 

Case  III-c  is  a  uniform  annular  beam.  Case  Ill-d  is  an 
annular ly-truncated  Gaussian  beam. 
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Pseudo- gray  scale  plots  of  transmitter  and  focal-plane 
irradiance  distributions  for  the  three  beams  in  Fig.  4.  A 
moderate  amount  of  thermal  blooming  is  present  in  the 
propagation  path. 
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FOCAL  PLANE  IRRADIANCE,  ARBITRARY  UNITS 


Figure  6.  Peak  focal-plane  irradiance  versus  total  transmitter  power 
for  the  four  beams  in  Figs,  5  and  6.  The  axis  units  are 
normalized  as  discussed  in  the  Appendix, 
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any  significant  advantages  for  reducing  thermal  blooming  when  the  goal  is 
to  achieve  maximum  focal-plane  irradiance  for  a  given  transmitter  powerp 
In  particular,  the  laser  modes  considered  here  that  produce  an  on-axis  null 
in  both  the  near-field  and  far-field  have  no  advantages  over  intensity  dis¬ 
tributions  that  have  the  null  only  in  the  near-field  of  the  transmitter. 
Significantly  less  blooming  does  occur,  however,  for  certain  annular,  cir¬ 
cularly  symmetric  laser  modes  and  for  initial  irradiance  profiles  that  are 
very  uniform  in  the  non- zero  regions  (a  uniform  annulus,  for  example). 

One  question  that  is  frequently  asked  is  whether  the  ^’holes’’  intro¬ 
duced  into  beams  by  unstable  resonators  should  be  eliminated  or  whether 
they  may  provide  some  advantage.  Although  we  have  not  considered  a 
sufficient  number  of  cases  to  make  a  definite  conclusion,  our  studies  to 
date  of  filled  beams  and  of  beams  with  a  central  hole  or  obscuration  indicate 
that  the  hole  makes  very  little  difference,  at  least  as  far  as  thermal  bloom¬ 
ing  is  concerned.  When  the  goal  is  to  minimize  thermal  blooming  dis¬ 
tortions  so  that  maximum  focal-plane  irradiance  is  achieved,  then  it  is  much 
more  profitable  to  force  the  initial  beam  intensity  profile  to  be  as  uniform  as 
possible  in  its  non-zero  regions  rather  than  attempt  to  produce  or  eliminate 
central  holes  in  the  beam.  Increasing  the  initial  transmitted  beam  diameter 

is  also  quite  effective  in  reducing  blooming  distortions,  particularly  when 

9 

coupled  with  adaptive  optical  techniques. 

Since  it  may  be  impractical  in  many  high  power  laser  applications  to 
achieve  either  sufficiently  large  apertures  or  to  produce  the  uniform  or 
special  intensity  distributions  that  minimize  thermal  blooming,  other 
irradiance-tailoring  techniques^^^  or  phasefront- tailoring  techniques  ^ 
may  be  more  useful.  Additional  work  is  now  in  progress  in  both  of  these 
areas. 
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APPENDIX 


In  this  appendix,  we  derive  the  normalizations  used  in  the  text  for 

« 

the  focal-plane  irradiance  and  for  the  field  ampjitudes  for  the  various  beams 
in  terms  of  the  total  laser  transmitter  power.  The  total  transmitted  power 
is  given  by 

'A-'' 

T 

where  is  the  aperture  area,  E  is  the  field  amplitude  as  defined  in 
Eq.  (1),  and  is  a  constant  given  by 


^l^^o^Po^  =  I  f^(^^2  p^x)  xdx  .  (A-2) 


The  maximum  aperture  dimension  is  and  is  the  e -folding  intensity 
radius  of  the  Gaussian  beam  in  Eq.  (  8). 

The  focal  plane  irradiance  as  plotted  in  Figs.  3  and  6  is  found 
from  Eqs.  (3)  and  (4): 


-  r 2tt  cos(me  )T 

I  U(r^,  e^,  z)  I  2 


(r  , 
m  o 


z) 


(A-3) 


with  I  defined  by  Eq.  (4).  We  are  interested  in  the  maximum  focal-plane 
m  ■' 

irradiance,  which  will  occur  at  6^  =  0  and  at  some  r^,  depending  on  the 
initial  irradiance  distribution.  Using  6^  =  0  combining  Eqs.  (A-2), 
(A-3),  and  Eq.  (4),  the  irradiance  can  be  written  as 


u  =  4'w(m  +  1 ) 


■  az 


FllV^o' 


(A-4) 


where  a  is  the  propagation  medium  absorption,  Fj  is  defined  in  Eq.  (A- 2), 


and  F2  is  given  by 


F2(a„/p„.  P)  = 


•a  / \^2  p 


f  (\^ 2  p  x)  J^(px)  xdx 
m  m 


{A-5) 


and 


2tt  n/Z 


P  = 


P  r 

o 


\z 


(A-6) 


We  have  evaluated  the  integral  function  in  Eq.  (A-5)  for  each  of 

the  7  radial  f^  functions  discussed  in  the  text  (Eqs.  (6)-(12))  to  find  the 

maximum  valTe  of  F^  as  a  function  of  the  normalized  parameter  p.  The 

1  /  -  1  d?  a  /b  -  7  r  /a  =  2/3,  and  p  =  1  were  used  for  this 

values  a^/p^  =  1.52,  a^/b^-  1,  ^rn'  o  . 

evaluation.  Table  A-1  gives  the  results  of  this  analysis.  To  obtain 

numerical  results,  we  chose  K=  10.  6  p.m,  z  =  2  km,  and  p^  -  23  cm 

(corresponding  to  a^  =  35  cm).  The  irradiance  plotted  as  the  ordinate 

axes  in  Figs.  (3)  and  (6)  is  thus  given  by 


u  I  ^  =  0.  0592  (m  +  1 )  (S.  R.  ) 


(A-7) 


in 


_  units  of  power  per  square  centimeter,  with  (S.  R.  )  being  the  strehl 
ratio  determined  by  the  thermal  blooming  computer  code  (S.  R.  =  1  for 
free-space  propagation).  For  the  blooming  calculations,  the  following 

parameters  were  taken: 


62 


absorption  =2x10  m 
transverse  wind  velocity  =  10  m/ sec 
beam  slew  rate  =  20  mrad/ sec 


Table  A-1.  Values  for  Peak  Focal- Plane  Irradiance 
Calculations 


Case 

m 

(F^)  max. 

(p)  max. 

Inf.  Gaus. 

0.  25 

0.  50 

0 

Trunc,  Gaus. 

0 

0,  225 

0.  342 

0 

Il-a 

1 

0.  225 

0.  167 

2.  50 

II -b 

1 

0.  508 

0.  279 

2.  31 

III- a 

0. 0509 

0.  161 

0 

ni-b 

0 

0.  508 

0.  536 

0 

m-c 

0 

0.  566 

0.  566 

0 

ni-d 

0 

0.  214 

0.  331 

0 
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